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The relevance of the Dirac equation for computations of nuclear structure is moti- 
vated and discussed. Quantitatively successful results for medium- and heavy-mass 
O nuclei are described, and modern ideas of effective field theory and density functional 

^ ; theory are used to justify them. 
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INTRODUCTION 
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To understand how to build atomic nuclei with the Dirac equation, we will begin by 
asking some simple questions. What are the basic nuclear properties that we are trying to 
correlate and predict? Why use hadrons (rather than quarks and gluons) as the degrees 
of freedom? Why use the Dirac equation rather than the Schrodinger equation to describe 
the dynamics? How can we build a simple model of nuclear matter that reproduces the 
empirical equilibrium properties and that can be extended to calculations of medium- and 
heavy-mass nuclei? How does the Dirac approach predict the nuclear shell model? And how 
can we relate the hadronic description of nuclei to the underlying strong-interaction theory 
of quantum chromodynamics (QCD)? 

The basic properties of nuclei provide stringent constraints on any nuclear theory. An 
accurate description of these properties is necessary for any useful predictions or extrapola- 
tions. We will concern ourselves primarily with bulk and single-particle nuclear properties, 
as listed below; a more detailed discussion can be found in Refs. Ill yl . 



X 

We certainly want to reproduce the observed shapes of nuclei: the interior density of 
a heavy nucleus should be relatively constant, since the nuclear forces "saturate" at the 
equilibrium density of nuclear matter (roughly p fa 0.15 fm -3 ). Moreover, the nucleus 
should have a well-defined surface, with the density decreasing from 90% to 10% of its 
central value over a distance of roughly 2 fm. Finally, because of saturation, the radius R of 
a nucleus should scale according to R ~ A 1 ^ 3 1.1 fm, where A = N + Z is the total number 
of neutrons (N) plus protons (Z). 

The total energy of the nucleus should follow the "liquid drop" formula 

E = - ai A + a 2 A 2/3 + a 3 Z 2 /A 1/3 + a 4 (N - Zf/A + ■ ■ ■ , (1) 



where typical values for the a, coefficients are given in Ref. [1|. 

The particle spectrum is determined by the qualitative features of the single-particle po- 
tential. In a nonrelativistic (Schrodinger) language, the central potential is midway between 
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a harmonic oscillator and a square well; this shape determines the ordering of the levels as 
a function of the orbital angular momentum. (See Ref. [![, Figs. 57.1 and 57.2.) In addi- 
tion, the spin-orbit potential is strong, which is instrumental in determining the major shell 
closures and, hence, the shell model. We will see below how these features are reproduced 
in a description based on the Dirac equation. 

These simple nuclear features are the ones we will focus on. We expect that they can be 
described adequately by a single-particle equation with an effective, one-body interaction. 
Such an approach has many names, depending on the system being studied and on the 
practitioner: "shell model", "mean-field theory", "Kohn-Sham" density functional theory, 
etc. Our goal is to correlate (fit) a modest number of nuclear bulk and single-particle data 
and then to predict other, similar data as well as possible. 



II. WHY USE HADRONS? 

Well, why not? Our focus is on long-range nuclear characteristics, and all measured 
observables are colorless. (In fact, most of the observables relevant to us are dominated by 
the isoscalar part of the interaction.) Moreover, hadronic variables (baryons and mesons) 
are efficient, since hadrons are the particles that are observed in experiments. Colored 
quarks and gluons participate only in intermediate states, and such "off-shell behavior" 
is unobservable; by using hadrons, we expend no theoretical effort combining quarks and 
gluons into color singlets that can actually be observed. 

So we pick the most efficient degrees of freedom by choosing hadrons. We will have 
to parametrize the nuclear hamiltonian anyway, since we cannot compute its true form 
from QCD, and hadronic variables, if combined in all forms consistent with the underlying 
symmetries, provide sufficient flexibility for our parametrization. We cannot guarantee that 
a single-particle hadronic approach will be successful in describing the observables of interest, 
but we want to see how well we can do. 



III. WHY USE THE DIRAC EQUATION? 

To motivate the Dirac equation as straightforwardly as possible, compare the particle 
spectrum (and fine structure) in a light atom with the spectrum in a heavy nucleus. An 
example of the former is given in Ref. 0], while an early example of the latter is given 
in Ref. jij, which is reproduced in Fig. 57.3 of Ref. pj. The most striking result is that 
it is impossible to draw the atomic fine structure to scale, since the splittings are roughly 
1/10,000 as large as the major-level splittings (at least for the deeply bound atomic levels). 
In contrast, the nuclear spectrum shows that the "fine" structure is really "gross"; the 
fine-structure splittings are as large as the major-level splittings to within a factor of two! 

The implication is that there must be some relativistic effects that are important in nuclei 
(unlike light atoms), and thus it is much more natural to use the Dirac equation to describe 
the quasi-particle nucleon wave functions. We will now try to understand this result by 
building a simple model of uniform nuclear matter. 



3 



IV. A SIMPLE MODEL OF NUCLEAR MATTER 



We consider a model first proposed by Walecka [5|, which contains nucleons (ip) and 
neutral (isoscalar) Lorentz scalar (0) and vector (V 1 ) mesons. This model is often referred 
to as "quantum hadrodynamics I" (or QHD-I, for short). The lagrangian density for this 
model (using the conventions of Ref. and suppressing counterterms for simplicity) is 



The included degrees of freedom are the minimal ones that will allow us to understand the 
qualitative features of the nuclear many-body system, which is our goal. We will describe the 
system in terms of Dirac quasi-particles moving in classical meson mean fields, an approxi- 
mation that we will elaborate on and justify later. Note that the baryon current (density) 
Tp^fixTp is conserved. 

It is important to emphasize that the Lorentz scalar and vector fields are effective fields 
that are introduced to parametrize the nucleon-nucleon (NN) interaction. The quanta of 
these fields never appear "on the mass shell" as real particles in any of the calculations dis- 
cussed here. They are analogous to the phonons that describe electron-electron interactions 
inside a metal. 

If one computes the NN interaction using one-boson exchange [purely for illustration, 
since the coupling constants g s and g v in Eq. (j2J) are large], one finds a short-range repulsion 
(from V 1 ) and a mid-range attraction (from <ft), which is characteristic of the NN force. 
Explicit pion exchange is of minor importance for our observables of interest; isoscalar, 
scalar and vector fields are dominant for the bulk and single-particle properties of heavy 
nuclei (and some multi-pion exchange is simulated by our effective fields anyway). 

Consider nuclear or neutron matter at zero temperature. We can treat the mesons at the 
mean-field level by taking 



where (V) = 0, since we assume that we are in the rest frame of the uniform matter. Note 
that O and V are constants. 

Why should mean meson fields yield a reasonable description of the system, since QHD-I 
is a strong-coupling theory? Our goal is to construct an approximate energy functional of 
the scalar (p s ) and baryon (p B ) densities and to fit the parameters in this functional to bulk 
nuclear properties. The mean meson fields give us a convenient way to do this, since they 
satisfy the mean-field equations: 



C = V(*7^ - M)^j + -(^0^0 - m s 2 2 ) - -(c^K - d u Vrf 



(2) 



(0> = 00 



(3) 



occ 




(4) 



occ 




(5) 




(6) 



M* 
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where "occ" signifies the occupied quasi-particle levels. (In infinite matter, we sum over 
states with both spin projections and with momentum k < k F , where k F is the Fermi 
momentum.) 

We solve these equations for stationary quasi-particle states; the problem is self- 
consistent, since p s = p s (M*) both determines and depends on the wave functions (H, 0]. 
The nuclear /neutron matter energy function(al) then becomes 

W = ^ Pi + ^§ (M - M*) 2 + A [ kF dt t 2 (t 2 + M* 2 ) 1 ^ , (7) 
where the baryon density is 



P* = 7d> (8) 



3tt 2 

and the isospin degeneracy is A = 2 for symmetric (iV = Z) nuclear matter and A = 1 for 
pure neutron matter [Z = 0). (Note that, by definition, the Coulomb force between protons 
is turned off.) 

One can now minimize the energy density £ with respect to p B to find the equilibrium 
point, and use the empirical equilibrium point of nuclear matter (density = p « 0.15 fm~ 3 , 
binding energy = e ~ 16MeV) to determine the two unknown ratios 

and , 

which are expressed more conventionally (and less dimensionally) as 



n 2 M 2 o 2 M 2 
C 2 9*M = 357 4 C 2 _ = 273 g (g) 



The resulting nuclear /neutron matter binding curves and the self-consistent effective mass 
M* as functions of the density are shown in Figs. 1 and 2 of Ref . 0| . The important features 
of these results are: 

• Symmetric nuclear matter is a self-bound liquid with an equilibrium point as defined 
above. This illustrates the "saturation" of nuclear forces. 

• Pure neutron matter is (generally) unbound at all densities. This reflects the positive 
symmetry-energy coefficient [0,4 in Eq. ((JJ] that enters when the number of neutrons 
and protons is different. 1 

• The nucleon effective mass at equilibrium density is roughly Mq m 0.6M. This shows 
that the scalar mean field is roughly —400 MeV at equilibrium; the corresponding 
vector mean field is roughly 300 MeV, and the two fields cancel to produce the relatively 
small binding energy of 16 MeV. We turn now to a discussion of this point. 



1 In QHD-I, this coefficient is too small. One must add a p meson, which couples to the difference of proton 
and neutron densities, to achieve an accurate result. See Ref. 
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What causes the nuclear matter saturation and the relatively small binding energy? Let's 
expand 8 / p B from Eq. (J7J) in powers of k f 
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(10) 



The lowest-order Lorentz scalar and vector contributions (which are proportional to p B ) set 
the scale for the large mean fields. [See Eqs. (jlj) and (J3J).] This scale is consistent with 
chiral QCD counting rules 0, H| , but these two terms cancel almost exactly in the binding 
energy, leading to an anomalously small remainder. However, they add constructively in the 
spin-orbit interaction, leading to appropriately large spin-orbit splittings in nuclei 0, 
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It is important to notice the different behavior of the vector and scalar interaction terms 
in Eq. (jlUJI . Whereas the vector interaction enters at only linear order in p B , the scalar 
interaction enters at all orders; moreover, the leading scalar term at every order in p B looks 
exactly the same, and they all add constructively. These terms are precisely what one 
gets by shifting the nucleon mass in the nonrelativistic kinetic energy term 3fcp/10M from 
M — > M* « M — g 2 p B /m 2 . These additional, repulsive, velocity-dependent interactions 
reduce the strength of the lowest-order, attractive scalar contribution and are crucial for 
establishing the location of the equilibrium point of nuclear matter. Thus the different 
behavior of the vector and scalar interactions leads to large relativistic interaction effects 
in the nuclear matter energy density. In contrast, the relativistic corrections to the kinetic 
energy (the nonleading terms in the first pair of square brackets) are indeed small; this is 
not where the important "relativity" is. 



V. MEAN-FIELD THEORY FOR NUCLEI 

We now compute the bulk and single-particle properties of atomic nuclei using essentially 
the same simple lagrangian discussed above. Our treatment follows that of Ref. [lpj, which 
is more than twenty years old, but which is still sufficient to illustrate the important points. 
We will discuss modifications and more modern treatments later. 

The basic idea is to allow the mean meson fields to be spatially dependent, and we 
will consider only spherically symmetric nuclei for simplicity. We again look for stationary 
quasi-nucleon states, and so the mean-field equations become p 

V 2 0o(r-) - m 2 s (f) (r) = -g s p s (r) , (11) 

VVo(r) - m 2 v Vo(r) = -g v p B (r) , (12) 

{-i a- V + g v V (r) + (3[M - g s <f> {r))}Tp a {x) = e Q ^ Q (x) . (13) 

These are coupled, nonlinear, differential equations that must be solved self- consistently. 
They are sometimes called Dirac-Hartree equations 0, U3] but are more accurately described 
as Kohn-Sham equations as we discuss in more detail below. 
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As one might expect, an accurate description of nuclear properties is not possible using 
only nucleons and isoscalar mesons. One must extend the model to include at least the 
Coulomb interaction between protons and an isovector p meson that allows for a more 
realistic description of the nuclear symmetry energy. (See Refs. [H Q for details.) The 
augmented model now contains four adjustable parameters: 

g s , g v , g p , m s . 

(The heavy meson masses are fixed at some "large" mass scale that is roughly equal to the 
nucleon mass M.) The couplings are fitted to the equilibrium point of symmetric nuclear 
matter and to the nuclear matter symmetry energy; the length scale, which is determined 
by m s , is set by fixing this parameter to reproduce the rms charge radius of a doubly magic 
nucleus, such as 40 Ca. 

Many nuclear structure calculations have been carried out within this relativistic mean- 
field theory (RMFT) framework. (See, for example, Refs. 0, or the extensive list of 



references in Ref. j6j.) One finds that the bulk properties of nuclei are well reproduced even 
in this relatively simple mean-field theory. Moreover, the single-particle spectrum reveals 
the well-known nuclear shell structure; this comes for free, since the parameters are fitted 
to the bulk properties of nuclear matter (and one nuclear length scale). 

Extensions of this simple model have been made to "fine-tune" the results. In the 1980's, 
numerous authors added terms involving nonlinear interactions of the scalar field: 

C= -^0 3 -^A0\ (14) 
and in the 1990's, various practitioners added vector self-couplings, like 

cr= +^c^ v 4 (^) 2 , (is) 

as well as other nonlinear and gradient-coupling terms, some motivated by the ideas of 
effective field theory; see the discussion below. (Many calculations in these extended models 
are cited in Ref. p. For an alternative approach that uses only nucleons in a lagrangian 
that contains numerous powers of fermion fields, see Ref. w]| and references therein.) 



These additional nonlinearities can be interpreted in terms of many-body nuclear forces, 
and they introduce additional density dependence into the nuclear energy functional, which 
allows it to more accurately reproduce the true energy functional. The new parameters 
are fitted either to additional nuclear matter properties, or to other theoretical calculations 
of nuclear matter (based on the Schrodinger equation), or to a selected set of data from 
finite nuclei. The basic conclusion from these extended calculations is that the successful 
qualitative features predicted by the original simple models persist, but the quantitative 
accuracy increased by nearly two orders of magnitude over a period of twenty years. For a 
comparison of the accuracy of results obtained with different collections of parameters, see, 
for example, Refs. 0,0]. For some recent state-of-the-art predictions of this approach (that 
is, calculations of nuclei that are not included in the fitting procedures), see Ref. (l5| . 

But it still remains for us to understand at a deeper level why these simple relativistic 
mean-field calculations can do such an excellent job of reproducing certain nuclear observ- 
ables. For this, we must study . . . 
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VI. MODERN DEVELOPMENTS 



The discussion in this section is a synopsis of the formalism presented in Refs. |6l 171 113111 
which is based on the ideas of modern effective field theory (EFT) and density functional 
theory (DFT). The interpretation of the earlier, successful results using EFT/DFT puts 
them on a firm theoretical basis. 

First of all, we interpret QHD as a nonrenormalizable EFT. This means that it contains 
known long-range interactions that are constrained by the underlying QCD symmetries, 
plus a complete (but non- redundant) set of generic short-range interactions, i.e., "contact" 
and "gradient" terms. The borderline between short and long ranges is characterized by the 
breakdown scale A of the EFT; empirically, we find that A « 600 M eV for QHD . 

If we ignore strangeness, then only nucleons and pions are "real" (stable) particles. The 
other field quanta are always virtual and just let us parametrize the NN interaction. As in 
any lagrangian theory, there are different ways to choose the generalized coordinates (fields) , 
but some coordinates may be more efficient than others jH,G3- 

The QHD EFT lagrangian explicitly exhibits the symmetries of QCD: The global, chiral 
SU(2)i x SU(2)r symmetry is nonlinear, approximate, and spontaneously broken p. The 
remaining global, isovector subgroup SU(2)y is realized linearly. It is straightforward (but 
usually tedious) to include electromagnetic interactions through the familiar local U(l) gauge 
symmetry 0]. 

The basic strategy for using the QHD lagrangian has been developed over the last several 
years 0, ; EB- First, assign an index v to each term in the lagrangian: 

v = d + n/2 + b , (16) 

where d is the number of derivatives (not counting those that act on nucleon fields), 2 n is 
the number of nucleon fields, and b is the number of non-Goldstone bosons. 

Now organize the lagrangian in powers of v and truncate. This gives an expansion in 
inverse powers of a heavy mass scale A « M, which has been shown to be reliable in 
calculations of medium- and heavy-mass nuclei [J El- Practically speaking, in the nuclear 
many-body problem, this expansion is in powers of k F /M, where k F is the Fermi momentum 
at equilibrium nuclear density (k F /M ps 1/3). 

Use the truncated lagrangian to construct an energy functional, which is to be interpreted 
within the DFT framework: We approximate the functional using factorized densities or 
fields, which produces a mean-field form of the functional. Expand it as a power series 
in density and momentum (by counting powers of v) and fit the remaining parameters 
to a restricted set of experimental data |8j]. These data typically include nuclear binding 
energies, prominent features of the nuclear electromagnetic charge form factors, and single- 
particle energy splittings for the least-bound orbitals 0. Define a set of Kohn-Sham (KS) 
single-particle orbitals that satisfy differential equations obtained by extremizing the energy 
functional with respect to the densities and fields. This procedure guarantees that all of the 
source terms in these equations are local. The KS orbitals are tailored to the generation of 
the ground-state density, and they include short-range and correlation effects adequate ly, i f 
the mean-field energy functional is a good approximation to the true energy functional 



[17]. 



2 Time derivatives acting on nucleon fields will generally bring down factors of the nucleon mass or energy, 
which are not small compared to A. 
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The mean-field energy functional constructed above omits some long-range contributions, 
which are generally nonlocal and nonanalytic functions of the densities. These contributions 
can be added systematically, by computing loop integrals using the well-known rules of EFT 
The effect of these loop contributions on the energy functional of atomic nuclei is an 
important topic for future study 

VII. SUMMARY 

The most important points in the preceding discussion can be summarized as follows: 

• The Dirac equation provides an economical and natural way to describe bulk nuclear 
properties and the nucleon single-particle spectrum, with the correct spin-orbit force 
(that is, the nuclear shell model) arising automatically. 

• Kinematical relativistic effects are small in nuclei, but dynamical relativistic effects 
from the interactions are important. 

• Modern QHD EFT's incorporate the basic symmetries of QCD. 

• The mean-field approach to heavy nuclei is really DFT, implemented through KS 
quasi-particle orbitals. The tested validity and accuracy of our truncation procedure 
for both fitted and predicted results shows that we really know something about the 
energy functional for cold nuclear matter near equilibrium density! 

• The energy functional can be extended beyond the mean-field parametrization using 
well-defined rules of EFT to compute the long-range contributions of loop integrals. 
This has been done recently [18l| . 

• The QHD/EFT/DFT/KS formalism provides a true representation of QCD in the 
low-energy nuclear domain. 

The basic message of this talk is: the Dirac equation is relevant for nuclear-structure 
physics, even though you might not expect it to be. Our quantitative successes justify its 
usage, but the modern theoretical ideas of EFT and DFT explain why it works. 
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